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Abstract. - We investigate one dimensional models (the Blonder, Tinkham, Klapwijk model 
and a tight-binding model) of non local transport at normal metal / superconductor (NS) double 
interfaces. We find a negative elastic cotunneling crossed conductance, strongly enhanced 
by additional scatterers away from the interfaces, suggesting the possibility of reflectionless 
tunneling non local transport at double NS interfaces with contacts having a sufficiently small 
extension. 



Introduction. - Single electron tunneling in a superconductor is prohibited if the applied 
bias voltage is smaller than the superconducting gap. However, an electron in the spin-up band 
can be reflected as a hole in the spin-down band, a phenomenon called Andreev reflection [1]. 
A charge 2e is transmitted in the superconductor at each Andreev reflection, so that the 
conductance of a highly transparent normal metal / superconductor (NS) contact is doubled 
compared to the one of the corresponding NN contact. The equilibrium properties of the 
superconductor (such as the value of the self-consistent superconducting gap) are modified 
by a normal electrode connected to the superconductor, a phenomenon called the inverse 
proximity effect. It is expected that most of the inverse proximity effect takes place on a 
length a if the area of the contact a 2 is much smaller than the superconducting coherence 
length £ [2] . The influence of the inverse proximity effect on transport properties can then be 
neglected, and a single channel, ballistic, one-dimensional model with a step- function variation 
of the superconducting gap captures the essential physics of localized interfaces, as shown by 
Blonder, Tinkham, and Klapwijk (BTK) [2]. Moreover, BTK introduce a repulsive potential 
at the NS interface, characterized by the dimensionless parameter Z, being the strength of 
the repulsive potential normalized to the Fermi energy. Transparent interfaces correspond to 
Z = and tunnel interfaces correspond to Z >• 1 . 

Disorder in the normal metal modifies strongly subgap transport at a single normal metal / 
insulator / superconductor (NIS) interface [3,4]. The conductance can be enhanced by orders 
of magnitude by constructive interferences in which an electron can "try" the tunneling process 
a huge number of times [3]. This effect due to scattering by disorder is already present in 
simple double barrier one-dimensional models. Melsen and Beenakker [5] consider a NINIS 
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Fig. 1 - Schematic representation of the electrical circuit corresponding to the NISIN junction studied 
experimentally in Ref. [22]. The current I a through electrode "a" is determined in response to a 
voltage Vb on electrode "b", with V a = 0. 



double junction in one dimension, with, in the BTK language, barrier parameters Z\ (for the 
NIN interface) and Z% (for the NIS interface). The conductance, averaged over the Fermi 
oscillations, shows a maximum for a value of Z\ comparable to Z 2 as Z\ increases while 
Z 2 ^> 1 is fixed [5]. This enhancement of the conductance shows that the double barrier 
model captures multiple scattering as in a disordered system. 

We address here similar effects for non local transport in NISIN junctions [6-20] in which 
the normal electrode "a" is at potential V a , the electrode "b" is at potential Vb, and the 
superconductor is at potential Vs (see Fig. 0. The non local conductance G a ,b(Vb) contains 
the information on how the current I a (Vb) in electrode "a" depends on the voltage Vb on 
electrode "b": Q a ,b(Vb) = dI a (Vb)/dVb [9,13]. The superconductor is taken as the reference 
voltage (Vs — 0), and we focus on the case V a = 0. Such devices have been realized in two 
recent experiments, performed in Karlsruhe by Beckmann et al. with ferromagnets [21], and 
in Delft by Russo et al. with a NISIN trilayer [22]. A sizeable crossed signal is measured in 
the latter [22] , which is surprising in view of lowest order perturbation theory in the tunnel 
amplitudes predicting an exact cancellation between the electron-electron and electron-hole 
channel crossed conductances [13]. We take as a working hypothesis that non local transport 
with normal metals is described by higher order contributions in perturbation theory in the 
tunnel amplitudes. These were already evaluated in Ref. [17] within microscopic Green's 
functions for localized interfaces. This approach was continued in Ref. [18] to account for 
extended interfaces with a large normal metal phase coherence length, giving rise to weak 
localization. Our task here is to investigate related issues in simple one dimensional models 
in the spirit of Ref. [5] . 

Blonder, Tinkham, Klapwijk (BTK) approach to a NISIN junction. - Let us first consider 
a one dimensional model of NISIN double interface within the BTK approach [16] (see Fig.Et). 
The gap of the superconductor is supposed to have a step-function variation: A(z) = A8(z — 
R/2)6(R/2 — z), and we suppose J-function scattering potentials at the interfaces: V(z) = 
H8(z + R/2) + H5(z — R/2) [2]. The two-component wave-functions are given by 

Mz) = ( I ) z lkFZ + a (i) elkFZ + b ( I ) e ~ lkFZ W 



S. Duhot and R. Melin: Reflectionless tunneling 



3 



N - R /2 S R/2 N 
(a) 1 1 



N -L/2 N -R/2 S R/2 N L/2 N 
(b) 1 1 1 1 

Z' Z Z Z' 



(c) 
(d) 



N -R/2 s R/2 N 



► z 



t t 

N -L/2 N -R/2 S R/2 N L/2 N 



t' t t t' 



Fig. 2 - (Color online) Schematic representations of the one dimensional models: the BTK model for 
NISIN (a) and NINISININ junctions (b), and the tight-binding model for NISIN (c) and NINISININ 
(d) junctions. 

UO \ ih F {z+R/2) p -{z+R/2)/(, , j / V \ -ik F {z+R/2) -(z+R/2)/£ 



UO \ e -ik F {z-R/2) e (z-R/2)/$ +d > ( V \ ik F (z-R/2) e (z-R/2)/i 

va J \ u o J 

ij 3 (z) = «'(?) e"^-^) +6' ( J ) e ^(»-H/2) } (3) 

where ipi(z), ip2(z) and ^(z) correspond respectively to z < —R/2, —R/2 < z < R/2 and 
R/2 < z, and Uq = 1 — v 2 = (l + i\/A 2 — w 2 /uS) /2 are the BCS coherence factors. We 
introduce the parameter Z — 2mH /fi 2 kp. The unknown coefficients a, b, a', 6', c, d, c', d' are 
determined from matching the wave- functions and their derivatives [2]. Assuming R £, we 
expand a' and b' to first order in exp (— to find the transmission coefficients 

jf" ^^|a'(fcW?)| 2 = (^i - ^6 + ^ + -) exp (- 2i? /^ + ° (° xp (- 4 ^)X 4 ) 
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.,_ \V(k F R)\ J = (7^1-^ + 7^ + -) exp(-2fl/0 + O(exp(-4R/0)(5) 



at w = 0. We deduce the first non vanishing term in the large- R, large- Z expansion of the 
non local transmission: 

T' = r^f^(\a'(k F R)\ 2 -\b'(k F R)\ 2 ) (6) 

JO ^7T 

= -^§8 exp (-2R/0 + O (exp (-4R/0) , (7) 

having a sign dominated by elastic cotunneling, in agreement with the Green's function ap- 
proach in which the first term in expansion of the non local transmission appears at order 
T 4 exp (— 2i?/£), where the large- Z normal transmission coefficient is proportional to Z~ 2 [2]. 

In the case of highly transparent interfaces corresponding to Z = 0, we find no crossed 
Andreev reflection: a'(u>) = 0, in agreement with the Green's function approach in Ref. [17]. 
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Fig. 3 - (Color online.) Variation of the crossed conductance G a ,b (in units of e 2 jh) for the junction 
on Pig.|5p, with Z' = Z\ = Za and Z = Z-z = Z^, — 10. (a) ... (f) correspond to an increasing values 
of the precision in the evaluation of the Fermi phase factors related to the superconductor. We have 
also shown the much smaller crossed conductance of the NINISIN junction, as a function of Z\ for 
the NIN contact, with the same value of Z for the NIS contacts. 



The elastic cotunneling transmission coefficient for Z = is given by 
1^12 _ 4(A 2 -^)exp(-2i?/g) 

1 { U 3w 2 [l-exp(-2i?/£)] + A 2 [l + exp(-2i?/£)/2]' 1 ' 



NINISININ junction. To describe multiple scattering in the normal electrodes, we 
consider now two additional scatterers at positions zi — —L\/2 in the left electrode and 
Z2 = L2/2 in the right electrode, described by the potentials V'(z) = H'8(z — zi) + H'S(z — Z2), 
and leading to the barrier parameter Z' = 2mH' /h 2 kF (see Fig. 03 for the definitions of Z and 
Z'). We average numerically the non local transmission coefficient over the Fermi oscillation 
phases <pi — kp{R — L{)/2, tp = kpR and (f2 = fc_p(i? — L 2 )- 

The variations of the crossed conductance at zero bias as a function of Z' for a fixed Z are 
shown on Fig. as well as the corresponding crossed conductance for the NINISIN junction. 
The integration over the microscopic Fermi oscillation phases for the latter involves a double 
integral so that the accuracy is larger than for the NINISININ junction involving a triple 
integral. As it is visible from the curves (a) - (f) on Fig. [31 corresponding to an increasing 
precision in the evaluation of the integrals, the crossed conductance for Z\ = has not 
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Fig. 4 - (Color online.) The notations used in the evaluation of the Green's functions of the one 
dimensional tight- binding model on the segments [ct,/3] (b), from the Green's functions of the full ID 
chain (a). 



converged to the limiting value obtained for the NINISIN junction, meaning that the change 
of sign in the crossed conductance at small Z\ for the NINISININ junction is an artifact related 
to the lack of precision in the evaluation of the triple integral (the crossed conductance at 
Z\ = for the NINISIN junction is indeed negative). The variation of the crossed conductance 
on Fig.OHshows a strong enhancement by the additional scatterers, suggestive of reflectionless 
tunneling, as for a NIS interface [5]. 

Green's functions. - Now we consider the same one dimensional geometry within Green's 
functions, and first evaluate the normal and superconducting Green's functions with appropri- 
ate boundary conditions. In one dimension, the Nambu Green's function of a superconductor 
at distance R and energy u) is given by 



g 1A (R,u) gi i2 (-R,w) 



(9) 



with 



g 1A (R,uj) = 
52,2 (R,v) = 



1 

2T 
1 

2T 



• cos (k F R) + sin (k F R) 



cos (k F R) — sin (k F R) 



g 1<2 {R,uj) = 52,i(-R,w) 



1 A 

2r7 



cos (kpR)e 



(10) 

(11) 
(12) 



with s = V A 2 — lu 2 and £,(lu) — Tlvf/s, where T is the bulk hopping amplitude of the one 
dimensional tight- binding model, and vf the Fermi velocity. 

The Green's functions on the finite segment [a, [3] can be deduced from Eq. by in- 
troducing a self-energy that disconnects the chain [23] . With the notations on Fig 0] we 
find 



1,1 
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Fig. 5 - (Color online.) Variation of the crossed conductance G a ,t (in units of e 2 /h) within the 
tight-binding model for the junction on Fig. |5Jl, as a function of t' for t/T = 0.0499 (corresponding 
to Z = 10 in the BTK model). The curves (a), (b) and (c) correspond to an increasing precision in 
the evaluation of the integral. 



with 

4 



1 + e - 4R /^ - 2 cos [2k F R)e- 2R/ ^ ) ] . (15) 



Similar expressions are obtained for \ and g a . a - 

Fig-Elshows the Green's functions result for the variation of the crossed conductance of the 
NINISININ junction as a function of t' for a fixed t (see Fig- EJi) ■ The numerical convergence 
is much faster than for the corresponding BTK model calculation because of the reduced 
dimension of the matrix to be inverted. We find the same feature as for the BTK model: 
the crossed conductance is enhanced by additional scatterers, as in reflectionless tunneling. 
Imposing the same normal conductance in the BTK and in the tight-binding models leads to 
the relation 

Z - 1 - { " T) \ (16) 
2t/T ' v ; 

leading to a good (but not perfect) agreement for the crossed conductance when the tight- 
binding and BTK results are rescaled on each other. 

Conclusions. - To conclude, we have investigated simple one-dimensional models consist- 
ing of NISIN double interfaces, with additional scatterers away from the two interfaces, in the 
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spirit of Ref . [5] . We find a strong enhancement of the crossed conductance by the additional 
scatterers, suggesting that non local transport at localized double NIS interfaces is enhanced 
by orders of magnitude, like in reflectionless tunneling at a single NIS interface. The geometry 
studied here is such that the Thouless energy associated to the dimension of the structure 
parallel to the interfaces is larger than the bias voltage. This reflectionless tunneling regime is 
not expected to correspond to the experiment in Ref. [22] because of the extended interfaces 
in this experiment, but may be probed in future experiments with disordered normal metals 
and interfaces of reduced extension. Finally, we also evaluated the crossed conductance as a 
function of energy, and found no sign change when the energy is increased for the BTK model: 
the crossed conductance is dominated by elastic cotunneling at all energies. 

* * * 

The authors thank D. Feinberg and M. Houzet for helpful discussions and in particular D. 
Feinberg for participating in the early BTK model calculations for the NISIN structure. 
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